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ABSTRACT 

Anisotropic distributions of charged particles including two-stream distributions give 
rise to generation of either stochastic electric fields (in the form of Langmuir waves, 
Buncman instability) or random quasi-static magnetic fields (Weibel and filamentation 
instabilities) or both. These two-stream instabilities are known to play a key role in 
collisionless shock formation, shock-shock interactions, and shock-induced electromag- 
netic emission. This paper applies the general non-perturbative stochastic theory of 
radiation to study electromagnetic emission produced by relativistic particles, which 
random walk in the stochastic electric fields of the Langmuir waves. This analysis takes 
into account the cumulative effect of uncorrelated Langmuir waves on the radiating 
particle trajectory giving rise to angular diffusion of the particle, which eventually 
modifies the corresponding radiation spectra. We demonstrate that the radiative pro- 
cess considered is probably relevant for emission produced in various kinds of astro- 
physical jets, in particular, prompt gamma-ray burst spectra, including X-ray excesses 
and prompt optical flashes. 

Key words: acceleration of particles - shock waves - turbulence - galaxies: jets - 
radiation mechanisms: non-thermal - magnetic fields 



1 INTRODUCTION 

Formation of collisionless astrophysical shocks, their interaction with ambient medium and each other are tightly coupled with 
two-stream instabilities giving rise to generation of either electric or magnetic fields or both in the shock vicinity. In the simplest 
kinetic version of the two-stream instability a tenuous electron beam excites Langmuir waves resonantly (i.e., at k = u> pc /vb, 
where k is the Langmuir wave vector, u} pc is the background plasma frequency, vt, is the velocity of the electron beam) on the 
linear stage of the instability. However, even though these small-scale Langmuir waves can initially be ge nerated by the wave 



parti cle resonance, non-linear wave- wave interactions will then transform the wave energy to larger scales (jKaplan fc Tsvtovich 



19731 ). Modern numerical models support t hat very small-scale random fields are generated initially at the shock front and 



then they evolve to larger and larger scales Uaroschek et al. 2004 . 20051 : Nishikawa et al. 20031 . 2005). Curiously, most of the 



work focuses on generation of stochastic magnetic fields due to Weibel or filamentation instabilities, although highly efficient 
charged particle acceleration occurring at the shock fronts necessarily requires generatio n of stochastic electric field as well 
( Bvkov fc Uvarovlll993l . ll999Ti . hdeed, numerical simulations of relativistic shocks (e.g.. [Nishikawa et al.|[2005t ) reveal very 



strong fluctuations of the electric charge p, which are tigh tly linked (kE = 4irp) with the longitudinal electric fluctuations. 
Quasilinear approach performed mainly analytically ( Bret et al. 20051 ) in 3-dimensional case shows that oblique modes 



dominate both purely longitudinal and purely transverse modes. This is true also for the case when a magnetic field parallel 
the beam velocity direction is included (jBret et al.ll2006l ) even though this magnetic field can completely suppress the purely 
transverse filamentation and Weibel modes. Both magnetic and electric fields are produced in the oblique modes and either 
electric or magnetic energy density can dominate depending on situation. Being excited these waves are subject of sophisticated 
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evolution due to interaction with each other and with charged particles 1 Silva 20061 ). These processes are deeply nonlinear 
and highly sensitive to the input parameters of the particular model. 

Accordingly, there can be regimes when magnetic (electric) inhomogeneities are dominant during t he entire ev o lution 
of the system and the effect of the electric (magnetic) inhomogeneities can be discarded. For example, iDieckmannl |2005) 
considered an ultrarelativistic regime of the two-stream instability when a broad spectrum of Langmuir waves (including 
non-resonant, large-scale, ones k <C Upe/vb) is excited and persist for a long time in the plasma. However, different regimes 
when initially dominant magne tic turbulen ce then gives a way to the electric turbulence and vice versa are also possible (for 
a more instructive review, see Silva 20061 ) . It would, therefore, be highly desirable to have observational tests allowing to 
distinguish between these different regimes, which can only come from adequate interpretation of electromagnetic emission 
recorded from an object. 

Surprisingly, the radiation arising as relativistic charged particles interact with stochastic electric fields of the Lang- 
muir turbulence has not been studied yet in sufficient detail, even though a number of i mportant and useful r esults have 
bee n obtained for the past 40 years. Initially, this radiative process was considered by iGailitis fc Tsvtovichl (| 19641 ) (see 
also Kaplan fc Tsvtovich 19731 ) for the case of small-scale Langmuir waves (k 3> cu pc /c, c is the speed of light) and then 
the theory was extended to t he case of the spatially uniform (k = 0) electric fluctuations with the plasma fr equency 



I Tsvtovich fc Chikhachev 



1969h. T he issu e of the characteristic frequency of this emission process was addressed bv iMelrosg 



Colgate! ( 1967 ) "electrostatic bremsstrahlung" in contrast to "magneto bremsstrahlung" , while 



i 19711). who call it after _ 

Schlickeiserl faooj) considered the total power of this proc e ss in a monochromatic approximation 



Getmantsev fc Tokarevl (119721 ) and lChiuderi fc Veltril (|1974l ) demonstrated that an ensemble of relativistic electrons in- 



teracting with the Langmuir turbulence would produce the emission with the same spectral index as for standard synchrotron 
radiation. The degree of polarization (for the case of highly anisotropi c, one-dimensional Langmuir turbulence) was found 
to be as large as for synchrotro n radiation in a uniform magnetic field I Tsvtovich fc Chikhachev 19691 : Kaplan fc Tsvtovichl 



19731 ; IWindsor fc Kellogg! 1974[). Therefore, one might conclude that these two emission processes are undistinguishable ob- 



servationally (e.g.. I Windsor fc Kelloggll 19741 ) . We believe, this conclusion is too straightforward. Indeed, for a narrow energy 
distribution of emitting electrons or for a broad (power-law) distribution with sufficiently sharp low- and/or high- energy 
cut-off, the specral shapes of the emissions can be remarkably different, enabling us to distinguish between various radiative 
processes. 

However, full description of possible spectral regimes of this radiative process has not been presented yet, although 
the general theoretical foun dation for this is readil y avail able (|Toptygin fc Fleishmanlll987l ; iToptygin et al.lll987l ; iFleishman 
2005, 2006b). In particular, IToptygin fc Fleishman! |l987l) clearly demonstrated that the emission in the presence of random 
electric fields be rries a lot of general sim ilarities to the emission in the random magnetic fields (so called diffusive synchrotron 
radiation, DSR, Fleishman! 2005 . 2006b ). although the exact expression for the electron scattering rate should be adjusted 
accordingly to properly take into account the polarization and dispersion of the Langmuir waves. The physics lying behind 
these similarities relates to the diffusive random walk of the emitting particle as it is being randomly scattered by either 
electric or magnetic irregularities. To emphasize that this diffusive motion is a key property to describe the emission correctly, 
we will refer to this emission process as "Diffusive Radiation in Langmuir turbulence" (or DRL) to be distinguished from its 
cousin DSR. So far, there is a number of terms suggested for this emission process, e.g., inverse Compton scattering of the 
plasma waves, inverse plasmon scattering, and electrostatic bremsstrahlung, although none of them is commonly accepted, 
since none of the titles is indicative enough for the radiative process considered. It is worth noting that the term "electrostatic 
bremsstrahlung" is not well suited for the emission in the Langmuir waves because a more usual "bremsstrahlung" is even 
more "electrostatic", than the DRL process discussed here. 



2 PERTURBATION THEORY OF DRL 



Perturbative treatment of electromagnetic emission by a charged particle assumes that the particle moves rectilinearly with 
constant velocity but takes into account non-zero acceleration of the particle in the external field. This perturbative treatment 
is widely used because of its simplicity. Typically, one calculates first the particle acceleration w(£) due to a given field along 
the rectilinear trajectory and then uses this expression obtained for w(t) to find the radiation spectrum. In the case of a 
random external field, however, when w(t) is also a random function of time t it is more convenient to express the radiation 
intensity via spatial and temporal spectrum of th e external e l ectric a nd /or magnetic field. 

Within theoretical formulation proposed by Fleishman! ( 2006b[ ) the spectral and angular distribution of the emission 
produced by a single particle with the Lorenz-factor 7 in a plasma with random field has the form: 



(2tt) 3 Q 2 
M 2 c 3 -y 2 V 



where 7* 



w 7; 



2u/ 2 7* 



dq dq5(uj' - q 



qv) 



(1) 



-1/2 



Q, M, and 7 > 1 are the charge, mass, and Lorenz-factor of the emitting particle, V is the 
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volume of the emission source, F„ 



is the Fourier component at the frequency go and the wave vector q of the Lorenz force 



transverse to the emitting particle velocity, 

, .2 



, UJ 
UJ = — 

2 



7 



(2) 



6 <C 1 is the angle between the wave vector of emitted wave k = kn and the particle velocity vector v, ui is the emitted 
frequency. Contribution W„ w (marked with the superscript _L) is provided by a component of the particle acceleration 
transverse to the particle velocity. In case of the electric (in contrast to magnetic) field, there is also a component of the 
acceleration along the particle velocity. The corresponding contribution has the form 

_ 2(2tt) 3 Q 4 



ir, 



dq dq8(uj' - q Q + qv) | E q 



(3) 



which is typically small by a factor 7~ 2 compared with the transverse contribution. Nevertheless, there exist special cases 
(e.g., a particle moving along one-dimensional turbulent electric field) when the transverse contribution is zero or very small 
and the parallel contribution comes to pl ay. For example, this is the ca se when a charge particle moves along the electric field 
of one- dimensional Langmuir turbulence IjFleishman fc Toptygin 2007). 

However, below we will consider only the transverse contribution, which is indeed the dominant one in most of the 
cases. Accordingly, the spectral distribution of the radiated energy is given by integration of {1} over the full solid angle 
dQ = sin 9d8dip « 2ird(uj' /u): 



(2^)V 
M 2 c 3 j 2 V 



d 



1/27? 



+ 



2w /2 7? 



q dq5(uj' - q + qv) | F, 



(4) 



similar to the DSR case |Fleishmanll2006bh . but with the Lorenz force F = QE specified by electric E in place of magnetic 
field. 



2.1 One- wave approximation in the DRL theory 



Let's consider first the simplest case when there is only one long Langmuir wave with fcpc <C ui ve , where kp is the w ave- vector 
and ojpe is the frequency of the Langmuir wave, which was described long ago by Tsvtovich fc Chikhachev ( 1969h . Here, the 
spatial and temporal spectrum of the electric field in this Langmuir wave takes a simple form 

TV 



E±(5d,q) \=A E 



(2tt) 



-S(q - k )(<5(go - w P c) + S(q Q +w pe )), 



(5) 



where Ae is specified by the energy density of the electric field in the wave. For example, for spatially uniform electric 
oscillations 

E = Eo COSWpet, (6) 

we have A E = E ± | 2 /4 . 

Calculation of the radiation intensity is extremely easy in case of field spectrum ([5}. Indeed, substituting ([5} into |l} 
and taking the integrals over frequency go and wave- vector q using the S- functions, we obtain the emission intensity (i.e., the 
energy emitted per unit time per unit frequency interval per unit solid angle) by dividing ([!]) over T: 



A E Q 4 
2nAPc 3 r 



1 + 



2^/2 



7 + ■ 



+ k v). 



(7) 



Let's discuss this emission intensity in more detail. Recall ([2]) that ui' is a function of the emission angle 0. On the other 
hand, because of the 5-function in (0, w' is a fixed number for a given set of values cu pc , kp, and v: 



kpv. 



(8) 



Therefore, there is strict correlation between the emission frequency and direction, thus, only one distinct frequency can be 
emitted along a given direction. 

Interestingly, emission intensity (J7J depends only weakly on the actual ko value as long as the condition fcpc <C lo pc 
holds. In particular, for a long Langmuir wave the radiation intensity is almost the same as for the spatially uniform temporal 
oscillations of the electric field with the plasma frequency. The radiation intensity into full solid angle (which is the differential 
spectral power, i.e., the energy emitted per unit time per unit frequency range) is given by simple integration of (JT} over 
d(u}' /w) with the use of the 5-function. It does depend on kp only weakly and for kp = has the form 
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in full agreement with the result of Tsvtovich fc Chikhachev (.1969). 

The dispersion of the plasma (terms Wp e /w 2 ) has almost no effect on the radiation intensity. The spectrum has unique 
asymptote 1^ cc lo 1 at low frequencies, lu <C UpcJ 2 with a peak at uj max = 2w pc 7 2 . The emission intensity vanishes at higher 
frequencies. 

Given that the radiation intensity for the case of single long Langmuir wave has only a weak dependence on the actual 
wavelength, it is tempting to extrapolate the presented results (in particular, spectrum (|9]) to the case when an ensemble of 
many long Langmuir waves exists in the source volume. Below we will see, however, that it is fundamentally incorrect, because 
the presence of the broad spatial spectrum of electric fluctuations makes a great difference compar ed with the one- wave case; 
general trends here are similar to those for the DSR in random magnetic fields (|Fleishmaidl2006bO . 



2.2 Perturbation theory for broad spatial spectrum 

Following the derivation given bv lFleishmanl (|2006ph . it is easy to find 
TV /„ v a v^ 



E„ 



(27T) 



K a p(qo,q) 



(10) 



where K a p{qo 1 q) = C a pK(q(,, q); C a /3 = q a qp/q 2 for the Langmuir turbulence since the electric field vector is directed along 
the wave vector in the Langmuir waves, while K(qo,<i) is the temporal and spatial spectrum of the Langmuir turbulence. 
Substituting (|10p into general expression ([TJ and dividing the result by T we obtain spectral and angular intensity of DRL 
per unit time, which is convenient to write down in the form: 



2tt 2 c 



(0 



+ 



2u/ 2 7? 



qL(oJ,9), 



(11) 



where qL(u>,8) = qn,{uj') is the effective scattering rate of the relativistic particle by the Langmuir turbulence, which plays a 
key role within the full non-perturbative treatment of the DRL: 



qL(uJ,9) = 



ttQ 2 



M 2 c 2 j 2 



a/3 



) K a 0(qo,<l)8(u' - q + qv). 



(12) 



To perform further calculations we have to specify the form of K a 0(qo, q). Initially, the streaming instability can give rise 
to highly anisotropic turbulent spectrum with a limiting case of purely one-dimensional turbulence, although randomization 
of the wave vectors will later result in more iso tropic turbulence patterns. D RL generated in the presence of one-dimensional 
Langmuir turbulence is specifically discussed bv lFleishman fe Topty gin (2007). By contrast, here we assume that the Langmuir 



wave vectors are isotropically distributed and the spectrum K(qo, q) can be approximated by a power-law above certain critical 
value ko: 

a»k" - 1 (El)q 2 

K ^o, q) = (fc g + g 2),/2+2 s (<*> - Wp.)- (13) 

Here, the presence of the <5-function is related to the assumption that the electric turbulence is composed of Langmuir 
waves all of which oscillate in time with the same frequency u pc ; the normalization constant a v is set up by the condition 
J K(qo, q)dgodq = (-E£), where (-EfA is the mean square of the electric field in the Langmuir turbulence. 

Now, substituting (|13[) into (|12|) . taking the integrals over dqo and dcos8 q with the use of two available (^-functions in 
(| 121) . and then integrating over dq, we find 




'■72 



+ k 2 



(14) 



where uj st = Q (El) 1/2 /Mc. Substitution of flT4) into (TTJ yields: 



2a u k^ 1 uj 2 t Q 2 
v{v + 2)c 2 7 2 



(0 



7 + 



LU 
2LJ 12 



7 +■ 



-v/2 



(15) 



Apparently, spectrum (|15|) looks rather differently from that in case of a single Langmuir wave (|9} . In particular, no 5-function 
enters (115[) . thus, a continuum spectrum rather than distinct frequencies is emitted along any direction. Clearly, there remains 
a distinct contribution to the emission intensity when lo' w lj pc . However, the range of the parameter space where this resonant 
condition holds is relatively narrow, so the "non-resonant" contribution from the remaining part of the parameter space where 
U)' 7^ Ljp C can easily dominate the resonant contribution. To see this explicitly, consider the radiation intensity into the full 
solid angle by integration of (115[) over dQ, that yields 
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Figure 1. Perturbative DRL spectra produced by a particle with 7 = 10 6 in a plasma with relatively weak (u> s t/k>pe = 3 ■ 10 — s ) long- and 
short- wave Langmuir turbulence for various LUo/uj pe ratios characterizing the turbulence: loq/u) pb = X, dashed/red curve; u)o/u) pe = 10 — 2 , 
dash- dotted /green curve; luq /u pe = 10~ 4 , d otted/bl ue curve. The spectrum of radiation in the presence of spatially uniform Langmuir 
oscillations tTsvtovich & C hikhachev\ \l96!\) is shown by solid/black curve for comparison. The spectra contain a distinct spectral peak 
at u) ~ 2o; pe 7 2 for long-wave turbulence, which is absent for the short-wave turbulence. 



2tt 



1/27, 2 



3ttc 



(16) 



where 



Ql(lu) 



o 2 7 V— 1 2 

37r a v k ui st 
2u(v + 2)c 3 7 2 7J 



1/272 



-2 . p 

7 H 



+ 



7 + 



i//2 



.(17) 



Here we express the radiation spectrum via the effective scattering rate = qL(u,d) averaged over the emission angle 

8, which defines the general non-perturbative expressions of the radiation intensity. Note that this averaging is performed 



with an appropriate weight described by the factor 



+ 



which enters Eqns (|ll|l and (|15[) . This factor is 



important since lu' is a function of the emission angle 9 according to equation ([2]). Combination of (|16|) and (|17[) yields finally: 

-1//2 



1 _ 47ra„fcp" 



d 



1/27? 



1 



7 + 
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2lo'' 



7 



+ 



(18) 



This integral cannot be taken analytically in a general case, but it is easy to plot corresponding spectra numerically. 
Figure [T] displays these spectra for a representative set of involved parameters. There are prominent differences in the DRL 
spectra in case of a broad spectrum of the Langmuir waves compared with the one- wave spectrum (oc to 1 ), which is plotted 
in the figure for comparison. Even though we cannot perform full analytical treatment of the spectrum, we can estimate 
the spectral shape in various frequency ranges. At low frequencies uj <C tj pc 7 2 , we can discard u)' in (I18|) everywhere in the 



braces except narrow region of parameters when uj' 



This means that for to UpcJ the integral is composed of two 



contributions. The first of them, a non-resonant one, arises from integration over the region, where u' <C u pe . Here, the 
emission is beamed within the characteristic emission angle of ■& ~ 7 -1 along the particle velocity. The integral converges 
rapidly, and so it may be taken along the infinite region, which produces a flat radiation spectrum, l w oc w°, or 1^ oc lu 2 at 
lower frequencies, u < ui pa "/. However, as far as u/ approaches oj pc , a resonant contribution comes into play. Now, in a narrow 
vicinity of u) pc , we can adopt 




-v/2 



oc 8(lo' 



(19) 



which results in a single- wave-like contribution, l u oc ui 1 . The full spectrum at lu < LOpef 2 , therefore, is just a sum of these 
two contributions. 

At high frequencies, w ^> Lu pc 'f 2 , the term lu' dominates in the braces, so other terms can be discarded. Thus, a power-law 
tail, Iuj oc lu~ v , typical for DSR high-frequency asymptote, arises in this spectral range, where there was no emission et all in 
case of a single Langmuir wave, «j2 . 1 1 
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Let's compare the DRL spectra with the DSR spectra in case of stochastic magnetic fields |Fleishmanll2006bl) . If the 
Langmuir turbulence consists of relatively small-scale waves, wo = koc ^ (j-> pa (see dashed curve in Figure [TJ then the shape 
of the spectrum is similar to the DSR spectrum. However, there is a remarkable difference in the case of the long-wave 
turbulence, coo <C w pc . Here, a distinct spectral peak at u) = 2w P c7 2 is formed with the linear decrease of the spectrum with 
frequency, there fore, the immediate vicini t y of th e spectral peak can indeed be described within the single- wave approximation 
as sugg ested bv lTsvtovich fc Chikhachevl |l969h . At lower frequencies, however, this falling part of the spectrum gives way to 
a flat spectrum, which is entirely missing within the one-wave approach. Position of the corresponding turning point depends 
on the u)o/u) pe ratio. It is worth emphasizing that the deviations of the DRL spectrum from the single- wave spectrum is 
prominent even for extremely long-wave turbulence, e.g., with loo/lo pc — 10~ 4 as in Figure [T] Therefore, the presence of a 
broad turbulence spectrum results in important qualitative change of the emission mechanism, which cannot generally be 
reduced to a simplified treatment relying on the single-wave approximation with some rms value of the Langmuir electric 
field. 



3 NON-PERTURBATIVE TREATMENT OF DRL 



The perturbative treatment of the emission produced by a relativistic particle moving in the presence of random fields breaks 
down sooner or later as soon as intensity of the field and/or particle energy increase (e.g., lFleishmanll2006bh . To find the 
applicability region of the perturbation theory applied above, we should estimate the characteristic deflection angle of the 
e mitting electron on the emission coherence length l c — 2cj 2 /lu, where the elementary emission pattern is formed. Similarly 
to lFleishmanl (|2006bl ) consider a simple source model consisting of uncorrelated cells with the size lo = 2-nc/uio, each of which 
contains coherent Langmuir oscillations with the plasma frequency ui pc . If ui p > u) ue then insid e each cell the electron velocity 



will change by the angle 60 ~ uist / (^07) , therefore, the consideration given in Fleishman] ( 2006b ) applies. However, if uiq <C uj pc 
then the electric field vector will change the direction approximately (uj pe /uJo) times during the time required for the particle 
to path through one cell, thus, the net deflection angle will be reduced by this factor (uj P!! /uiq) due to temporal oscillations of 
the electric field in the Langmuir waves, therefore 9o ~ w s t/(w P c7). Then, after traversing TV = l c /lo cells, the mean square of 
the deflection angle is 2 = 8qN ~ u^wo/(wWp e ). The perturbation theory is only applicable if this diffusive deflection angle is 
smaller than the relativistic beaming angle, 7 -1 , i.e., it is always valid at sufficiently high frequencies a> > w* = u 2 t uoj 2 /w pe . 
Note, that the bounding frequency u, increases with Wo, while DSR displays the opposite trend. Th e perturbation the ory will 
be applicable to the entire DRL spectrum if the condition 2 < 7 -2 holds for the frequency uj pe -y ([Fleishman! 12006 bf) . where 
the coherence length of the emission has a maximum. This happens for the particles whose Lorenz-factors obey the inequality 

(20) 



7 < CJ pe /Vs t U;o). 



We see, therefore, that generally, especially for relatively strong Langmuir turbu lence, the perturbative treatment is in- 
suffic ient to fully describe the radiation spectrum, thus, the n on-perturbative version dToptygin fc Fleishman 1987 ; Fleishman 
2005) of the theory should be explored. As demonstrated in iToptygin fc Fleishman! ()l987h the same general expressions for 



the radiation intensity produced in the presence of stochastic electric fields are valid like in the random magnetic fields, 
although the electron scattering rate by Langmuir waves qh{oj) = qL{u,9), which has already been introduced in the previous 
section, will differ from that in case of magnetic turbulence. As we will see below, all the differences between DRL spectrum in 
Langmuir turbulence and the DSR spectrum in magnetic turbulence are ultimately related to the difference in the expressions 
for the scattering rate qh[uS) for these two cases. 

Let's consider first the reg ime when there is no regular magnetic field and, thus, the radiation spectrum is defined as 
( Fleishman fc Bietenholzll2007n 
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where $(s) is the Migdal function (|Migdallll95l . Il956l ) 
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Ul 2 J 



(21) 



(22) 



dt exp(— 2st) sin(2sf) x 
having rather simple asymptotes for large or small s values: 

$(s) ~ 1, if s > 1, $(s) ~ 6s, if s < 1. (23) 

Figure [2] presents the DRL spectra for the case of relatively weak electric field, u> s t w pc ; the dash-dotted curves are the 
corresponding perturbative spectra. The non-perturbative effect (multiple scattering of the radiating electron by Langmuir 
waves) modifies the spectrum around the frequency aj pc 7 giving rise to asymptote I u tx uj 1 ^ 2 in this spectral region. Note that 
the frequency u>*, where the break from the oc a>° to oc u) 1 ^ 2 asymptote occurs, increases with ojo increase, while the DSR in 
the static random mag netic fields (|Fleishmanll2006bl ) displays the opposite trend. 
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Figure 2. DRL spectra produced by highly relativistic particle with 7 = 10 moving in a plasma with long-wave Langmuir turbulence. 
Full non-perturbative spectra are shown by solid curves, while the corresponding perturbative spectra are by dashed and dash-dotted 
curves. 

Figure 3. DRL in the presence of strong long-wave Langmuir turbulence. Dotted curves are calculated within the perturbation theory. 
Extremely strong suppression of the DRL spectra compared with the perturbative ones is evident. 



However, as shown in lSilva (|2006h the electrostatic field in the Langmuir waves generated at the shock front can be rather 
strong, e.g., of the order of nonrelativistic wave breaking limit, cu s t ~ w pc . In this case, the non-perturbative treatment is 
important at the full frequency range below the spectral peak at 2ui pc 'y 2 , Figure [3] For completeness of the possible DRL 
regimes considered, Figure [3] presents also the DRL spectra for the (less realistic) case of a very strong random electric field, 
ujst 3> Lo P c- Here, the non-perturbative spectrum deviates from the perturbative one even at the frequencies above 2lo pc j 2 , 



giving rise to a suppressed spectrum lu. 



->-l)/2 



a very broad non-perturbative power-law region of the spectrum, oc ui 1 ^ 2 , is formed. 



(compared with the perturbative one 1^ oc ui "). At the lower frequencies, 



4 DRL VS SYNCHROTRON RADIATION 



If both Langmuir turbulence and regular magnetic field are present in the source, then both DRL and synchrotron emission 
are generated. One could assume that the DRL is only relevant in ex treme conditions when the Langmuir turbulence energy 
density exceeds that of the regular magnetic field (e.g., Melrose! 1971 ). It is not the case, however, because these two emission 



mechanisms are efficient in differing frequency domains. 

Consider joint ef fect of the Langm uir turbulence and a regular magnetic field on the radiation spectrum. In this case full 



equations (30) from ( Fleishman 20051 ) together with the exact expression for 5l(i^) given by Eq. (|17[) should be used. It is 



worth emphasizing that the condition lu-r c <^ u> pc , where uJBe is cyclotron frequency in the regular magnetic field, is sufficient 
for full applicability of the available stochastic theory of radiation for the case under study. This condition means that the 
distortion of the electron trajectory due to the regular magnetic field is small during the period of electric oscillations, thus, 
it is not needed to be small at any spatial scale, so no further restriction on values of coo, ui st , and cu pc is necessary. 

Figure 2] shows examples of the spectra for the cases of short-wave or long-wave Langmuir turbulence superimposed 
on a regular magnetic field. The full spectrum consists of the standard synchrotron contribution (region 1^ oc uj 1 ^ 3 with 
exponential cut-offs at low- and high-frequency edges) and DRL contribution, which is the most prominent at the high 
frequencies lu 2> lobcJ 2 , although it is also present at sufficiently low frequencies, where the spectrum 1^ oc lu 2 is formed. 

Figure U allows direct comparison of the DRL spectrum with that of the standard synchrotron radiation. We point out 
that these two emission processes occupy distinct frequency ranges if ojb c "C ^ pc - In particular, the synchrotron spectrum 
displays exponential cut-off at the frequencies uj > c^Be7 2 , while the DRL displays here flat of even rising spectrum up to 
lu ~ LUp C "f 2 . This means that DRL can dominate this spectral range even if the energy density in the Langmuir turbulence is 
lower than the magnetic field energy density. This holds also for a power-law energy spectrum with a high-energy cut-off at 
some 7 m ax for the spectral range uj > 0JBe7max- 

Moreover, in sufficiently dense plasmas DRL can dominate the entire radiation spectrum even for the conditions when the 
Langmuir turbulence energy density is much smaller the the magnetic energy density. This happens for relatively low-energy 
(although ultrarelativistic) electrons, whose synchrotron emission is significantly suppressed by the Razin-effect (referred also 
to as density effect). Figure [5] displays such an example. Here the turbulent energy density is ten times lower than the magnetic 
energy density. Radiation spectrum produced by higher energy electrons (7 > 120 for the parameters selected to plot the 
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Figure 4. DRL spectra in the presence of uniform magnetic field with the energy density of the order of 10% of the Langmuir turbulence 
energy density. DRL dominates at frequencies above cjg e 7 2 , where the perturbative treatment is applicable as long as u s t *C tj pe . 

Figure 5. DRL spectra in the presence of uniform magnetic field. The Langmuir turbulence energy density is of the order of 10% of 
the magnetic energy density. Synchrorton spectrum is strongly suppressed by the effect of density. DRL dominates the entire spectrum 
for sufficiently low-energy electrons with ■y < 120. 



figure) consists of both DRL and synchrotron contribution, although the later becomes narrower and weaker as 7 decreases. 
Eventually, for 7 < 120, the synchrotron contribution disappears being strongly suppressed by the Razin effect in contrast to 
the DRL spectrum, which is less sensitive to the density effect. 



5 DISCUSSION 



Shock waves, in particular, jet shocks are believed to be the sites where various kinds of the two-stream instability can be 
operational. Depending on the conditions at the shock front and its vicinity either magnetic or electric fluctuations or both 
will be excited. Here, we specifically considered the case when random electric fields in a form of Langmuir wave turbulence 
dominate. 

Modern computer simulations of shock wave interactions, especially in the relativistic case, suggest that the energy density 
of the excited Langmuir turbulence can be far in excess of the energy of the initial regular magnetic field. In particular, at 
the wave front the electric field can be as strong as the corresponding wave-breaking limit, i.e., cu s t ~ uo pc . In this case the 
random walk of relativistic electrons in the stochastic electric field can give rise to powerful contribution in the nonthermal 
emission of an astrophysical object, entirely dominating full radiation spectrum or some broad part of it. 

So far, the D RL (electrostatic bremsstrahlung) has been applied to a number of astrophysical objects. For example, 
Schlickeiser noted that electrostatic bremsstrahlung is an attractive alternative to s tandard synchrotron r adiation to 

produce the observed nonthermal emission from jets in active galactic nuclei. In addition, ISchroder et al. (2005) developed 
a simplified model of the galactic diffuse sub-MeV emission based on monochromatic approximation of both synchrotron 
radiation and the DRL, which gives rise to a remarkably good agreement between the model and the observations. We point 
out that the use of presented here DRL spectra will be helpful to further develop that model especially in the range of 
high-energy cut-off of the radiation spectra. 

Although any detailed application of the considered emission process is beyond the scope of this paper, we mention 
that the DRL is also a promising mechanism for the gamma-ray bursts and extragalactic jets. In particular, some of the 
prompt gamma-ray bursts display rather hard low-energy spectra with the photon spectral index a up to 0. The DRL spectral 
asymptote I u oc lo 1 , which appears just below the spectral peak at 2io pe 'y 2 , fits well to those spectra. Remarkably, the flat 
lower - frequency asymptote, I w (x lo°, can account for the phenomenon of the X-ray excess jPreece et al.lll996l ; ISakamoto et al 
2005) and prompt optical flashes accompanying some GRBs. 



In addition, this mechanism (along with the DSR in random magnetic fields, iFleishm an 2006a) can be relevant to the 



UV- X-ray flattenings observed in some extragalactic jets. For example , although full spatially resolved radio to X-ray spectra 
of the jet in M87 agrees well with the DSR model ([Fleishman! l2006al ) . for the jet in 3 C 273 this agreem ent holds from the 



radio to UV range, while its X-ray emission seems to require an additional component ( Jester et al. 2006[ ). Alternatively, the 
entire UV-to-X-ray spectrum of 3C 273 might be produced by DRL, which can be much flatter than usual DSR (see, e.g., 
Figure [4| in the range WBe7 2 <C w <C o-> pc 7 2 . 
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